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CN ■ When several sparticle masses are known, the kinematics of SUSY decay pro- 

CJ ! cesses observed at the LHC can be solved if the cascade decays contain suffi- 

cient steps. We demonstrate four examples of this full reconstruction technique 
applied to channels involving leptons, namely a) gluino mass determination, b) 
sbottom mass determination, c) LSP momentum reconstruction, and d) heavy 
higgs mass determination. 
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Abstract 



1. INTRODUCTION 



The potential of the LHC for SUSY parameter determination has been studied in great detail for the past 
■ seven years Q|2l|3l|4lE|6]]. One of the most promising methods involves the selection of events from 

a single decay chain near the kinematic endpoint. Information on the masses involved in the cascade 
decay can be extracted from the endpoint measurements. It has been established that one can achieve a 
few percent accuracy for sparticle mass reconstruction using this technique with sufficient statistics. 
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In this paper we propose a new method for reconstructing SUSY events which does not rely only 
on events near the endpoint. Instead one kinematically solves for the neutralino momenta and masses of 
D . heavier sparticles using measured jet and lepton momenta and a few mass inputs. 



To illustrate the idea we take the following cascade decay chain 



X ■ 9 -» bb -> x U 2 bb - ebb£ - Xibbll (1) 

h ; 

This decay chain is approximately free from SM background with appropriate cuts. The five SUSY 
particles which are involved in the cascade decay have five mass shell conditions; 

2 2 

Al Al 

m2 X° 2 = ( P X? +Pfe) 2 > 

2 - (p x » +Ph + Pe 2 +Pb 1 ) 2 



nir 

o 



raj = (p^o +p h +p l2 +p bl +p b2 ) 2 - (2) 

Of these five masses, m^o ,m^o and can be measured at the LHC using first generation squark cascade 
decays with an accuracy of ~ 10% (the mass difference is measured more precisely). Moreover, with 
input from a future high energy Linear Collider these masses might be determined with an accuracy 
~ 0(1%). We therefore assume for the present work that the masses of the two lighter neutralinos and 
of the right handed slepton are known, and we ignore the corresponding errors. 

For a bbll event, the equations contain six unknowns (nig, m~ b and p^o) which satisfy five equa- 
tions. For two bbll events, we have ten equations while we only have ten unknowns (two neutralino 
four momenta, m- g and m^o ). Mathematically, one can obtain the sbottom and gluino masses and all 
neutralino momenta if there are more than two bbll events. 
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Table 1: Some sparticle masses in GeV at SPSla. 



We call this technique the "mass relation method" as one uses the fact that sparticle masses are 
common for events which go though the same cascade decay chain. Note events need not be near the 
endpoint of the decay distribution to be relevant to the mass determination. In the next section we 
demonstrate the practical application of this method to measurement of the masses of the gluino and 
sbottom. 

As a byproduct of the technique, once the mass of the squark and of all the sparticles involved in 
the decay are known, the momentum of the lighter neutralino can be fully reconstructed, and this further 
constrains the event. 

In SUSY events sparticles are always pair produced and there are two lightest neutralinos in the 
event. If squark decays via q — » %2 ~~ > & ~ > X? can be identified on one side of the event then the 
neutralino momentum can be reconstructed as described above. The transverse momentum of the lightest 
neutralino in the other cascade decay can then be obtained using the following equation 

Pt(x?(2)) = PT(miss) + p T (Xi(l)), (3) 

provided that there are no hard neutrinos involved in the decay. This transverse momentum can be used 
to constrain the cascade decay of the other sparticle. 

For the case where the other squark decays via q — * xtl ~ > XiQW followed by W — > q'q", the 
chargino mass can be determined by using Eq. © and the following relations, 

Pg = Px'l(2)+Pj + Pw, 

P% = m§, (4) 

where pj is the momentum of the selected high p? jet which comes from the squark decay and pw 
is the momentum of the two jet system consistent with the W interpretation. The neutralino momen- 
tum resolution is important for the chargino mass reconstruction and we discuss this in section 3. The 
reconstruction will be discussed in a separate contribution! 7 ]. 

The full reconstruction technique can be extended for higgs mass reconstruction. In section 4, we 
discuss the heavy higgs mass determination from the process H — ► X2X2 followed by jify—tll—* Uxi- 
This process is also useful for discovery of heavy higgs bosons. The four lepton momenta and missing 
momentum can be used to reconstruct the higgs mass assuming that the px of the higgs boson is very 
small. 



2. GLUINO CASCADE DECAY 

We first discuss the results of a simulation study of the process where a gluino cascade decays into a 
sbottom at model point SPSla[8|. The relevant sparticle masses for this study are listed in Table 1. The 
events were generated using the HERWIG 6.4 generator [9] [ 10 1 and passed through ATLFAST [11 ], a 
parametrised simulation of the ATLAS detector. 

We study only events which contain the cascade decay shown in Eq.®. We then apply the fol- 
lowing preselections to reduce backgrounds: 

• pjp iss > 100 GeV 

• M eff > 600 GeV 

• at least 3 jets with p T1 > 150 GeV, p T2 > 100 GeV and p T s > 50 GeV 



• exactly two jets with pt > 50 GeV tagged as 6-jets 

• exactly two OS-SF leptons with pm > 20 GeV, pxn > 10 GeV, and invariant mass 40GeV< 
ma < 78 GeV 

The solution of Eq. © can be written in the following form: 

mj = Fo + F im l±F 2 D, 
where D 2 = D + D x m\ + B 2 m\. (5) 

Here Fi and Di depend upon pg. and and the neutralino and slepton masses. In the event, there are 
two b jets and we assume that the b jet with larger pt originates from the b decay. The two leptons must 
come from \ 2 an ^ ^ decay. There are maximally four sets of gluino and sbottom mass solutions together 
with two lepton assignments for each decay, because we cannot determine from which decay the lepton 
originates. To reduce combinatorics we take the event pair which satisfies the following conditions: 

• Only one lepton assignment has a solution to the Eq. ® 

• For a pair of events there are only two solutions and there is a difference of more than 100 GeV 
between the two gluino mass solutions. 




Fig. 1: nig obtained by using Eq. J5} for two bbii events. 

In Fig 1 , we plot the minimum m g solution which satisfies the conditions given above. The peak 
position is consistent with the gluino mass, and the error on the peak position obtained by a Gaussian fit 
is around 1.7 GeV for 100 fb -1 . For the events used in the reconstruction, each event is used on average 
five times. Note that the a of the Gaussian fit is large (~ 56.7 GeV) and is determined by the resolution 
on the momentum measurement of the four 6-jets. It is worth stressing that the results presented here 
were produced by using a parametrised simulation of the response of the ATLAS detector to jets, based 
on the results of a detailed simulation. Results which crucially depend on the detailed features of the 
detector response, such as the possibility of discriminating the two sbottom squarks (see below) need 
to be validated by an explicit detailed simulation of the detector performed on the physics channel of 
interest. We only attempt here to evaluate the impact of the new technique on sparticle reconstruction. 

Once the gluino mass has been determined one can reconstruct the sbottom mass by fixing the 
gluino mass to the measured value. Here one needs to solve only Eq.©, which involves only two 6-jets 
in the fit, and therefore errors due to the jet resolution are expected to be less than those for the gluino 
mass reconstruction. 

For each event, there are two sbottom mass solutions m^(soll) and mr(sol2), each sensitive to 
the gluino mass input. The difference between the gluino and sbottom mass solutions is however stable 
against variation in the assumed gluino mass. The mass itself may have a large error in the absolute 
scale, but the mass differences are obtained rather precisely, as is the case in the endpoint method. 



In Fig. 2 (left), we plot the solutions for all possible lepton combinations in the m g — m~ h (soil) 
m g — m~ b (sol2) plane. Here we use the 6-parton momentum obtained from generator information. One 
of the solutions tends to be consistent with the input sbottom mass. Moreover the two decay modes 
g —> b\b and 626 are clearly separated. 

We can compare the results from the previous analysis with those from the endpoint analysis! 12 1, 
where one uses approximate the formula 

Pv«= I— P«- (6) 
* 2 V ma J 

This formula is correct only at the endpoint of the three body decay X2 ~~ > Xi^> but is nevertheless 
approximately correct near the edge of \2 ~ > ^ ~ > ^Xi f° r SPSla. The sbottom mass obtained by 
using Eq.© is shown in Fig. 2(right). For this case, the 62 peak at 70.6 GeV is not separated from the 61 
peak at 103 GeV. 




m(gluino)-m(sbottom,sol1) mgl-msb(GeV) endpoint approximation(parton level) 




Fig. 3: As for Fig. 2 but with the b jet momentum used instead of the b parton momentum. 



The b\ mass, or the weighted average of the sbottom masses, is easily obtained. The b jet resolution 
is not sufficient however to clearly separate the b\ and b\- This can be seen in Fig. 3 where the plots show 
the distributions corresponding to Fig.2(left) and (right) but now with the b parton momenta replaced by 
b jet momenta. For the endpoint analysis (Fig. 3 right), a correct evaluation of the sbottom masses would 
require a fit taking into account the shape of the response of ATLAS to b-jets. In order to approximately 
evaluate the achievable statistical precision, a naive double gaussian fit was performed on the distribution 
shown in Fig. 3 right, which corresponds to J dtL = 300 fb _1 . The resulting statistical uncertainties are 
±1 GeV (±2.5 GeV) for the m- g — m- b Jjrig — ) peak positions respectively. Additional systematic 
uncertainties, not yet evaluated, as well a 1% error due to the uncertainty on the jet energy scale should 
also be considered. These numbers are obtained assuming the presence of two gaussian peaks in the data. 

For the mass relation method the number of events available for the study is larger by a factor of 2 
because events away from the endpoints can be used. We also use the exact formula for the mass relation 
method. Although the analysis is more complicated due to the multiple solutions, we believe it to be a 
worthwhile technique for use when attempting to reconstruct the b\ and 62 masses. 

3. NEUTRALINO MOMENTUM RECONSTRUCTION 

In this section, we discuss the reconstruction of the momentum of the lightest neutralino. As we have 
discussed already, the mass shell condition can be solved for long decay cascades, such as q — > X2Q ~ > 
£q£ — > XiQ^- F° r this process we have two neutralino momentum solutions for each lepton assign- 
ment. One may wonder if the solutions for the neutralino momentum might be smeared significantly, 
because of the worse jet energy resolution as compared to leptons, and the jet px is generally much 
larger than the neutralino momentum for the cascade decay. In Fig. 4(left) we show the distribution of 
PT(reco)/pT (truth) for the point studied in Q. Here we choose the correct lepton combination using 
generator information, and take the solution which minimizes |pr(reco)/p^ (truth)— 1|. Except for the 
case where we took the wrong jet as input the reconstructed pt is within 20% of the true neutralino 
momentum. The result for the gluino cascade decay into sbottom Eq.(l) is similar. 

In Fig. 4(right) we show a similar reconstruction for the gluino cascade decay, but unlike Fig.4(left), 
we use both lepton combinations. We fix the gluino mass to the input value 1 and take events where one of 
the four sbottom mass solutions is consistent with the input sbottom mass such that — m^(best)| < 
10 GeV. We then take the solution where the sbottom mass is closest to the input mg . There are still two 
p^o solutions, and we choose the one which minimize | min(pT (reco)/pT (truth), pt (truth) /pr(reco)) — 
1 1 . The neutralino momentum resolution is worse than that obtained using the correct lepton assignments 
only. Nevertheless a significant fraction of events are reconstructed with 0.8 < |p(reco)/p(truth)| < 
1.2. 

4. HIGGS MASS RECONSTRUCTION 

A promising decay for the observation of heavy and pseudo-scalar higgs bosons in the difficult region 
with intermediate tan j3 is the decay into two neutralinos. When both neutralinos decay through the 
chain 

the resulting signature consists of events with four isolated leptons (paired in opposite-sign same-flavour 
pairs) and no jet activity. The main SM backgrounds to this signature are ti production, where both 
the 6-jets and the Ws decay into leptons and Zbb production. The key element for the rejection of 
these backgrounds is the fact that the leptons from b decays are not isolated. A detailed study of the 
performance of lepton isolation in the detector is needed to assess the visibility of the signal. Additionally 
there is an important SUSY background, including irreducible backgrounds from direct slepton and 

'Here we adopt an event selection which makes use of the true (input) gluino and sbottom mass values, although in practice 
fitted values would be used. 




p T (reco)/p T (truth) p|(reco)/|p| (truth) for the best solution 



Fig. 4: Left: The calculated \\ transverse momentum divided by the true transverse momentum. The decay — > X2Q ~ * 
tqt — » xiqtt is studied for the model point used for the chargino study: mo = 100 GeV, m 1( / 2 = 300 GeV, Ao = —300 GeV, 
tan/3 = 6, and /i > 0. Only the correct lepton choice is used. Right: |p(reco)/p(truth)| for the decay chain Eq.Q for SPSla. 

gaugino decay. Full background studies as a function of the SUSY parameters were performed by the 
ATLAS and CMS Collaborations fl3l [l4l . We propose here, along the lines of the previous sections, 
a technique for the complete reconstruction of the higgs peak, based on the knowledge of the masses 
of X2, £r and \i- I n this case one has 8 unknown quantities: the 4-momenta of the two LSP's, and 8 
constraints: six on-shell mass constraints (3 for each leg), and the two EJp iss components. 

To demonstrate the power of the method, we apply it to Point SPS la, for which the mass of the A 
and of the H is ~ 394 GeV. The BR into X2X2 is 6% for the A (H). We perform the study on 1000 
events for 

A - X2X2 - UUXixl 

corresponding approximately to the expected statistics for 300 fb _1 . We simply require 2 isolated leptons 
with px > 20 GeV and 2 further isolated leptons with p? > 10 GeV, all within \rj\ < 2.5. The efficiency 
of these cuts is ~ 60%. 

We have not performed any background simulations because at this stage we only wish to explore 
the viability of the full reconstruction technique. The main problem for the reconstruction is the correct 
assignment of the leptons to the appropriate decay chain. The first selection is based on requiring a 
unique identification of the lepton pairs coming from the decays of the two %2 S - We therefore require 
that either of the following two criteria is satisfied: 

• the flavour configuration of the leptons is e + e + r u + /x~ 

• the lepton configuration is either e + e~e + e - or / u + fi + fj," , but for one of the two possible pair- 
ings the invariant mass of one of the pairs is larger than 78 GeV, i.e. above the lepton-lepton edge 
for the X2 decay. 

The total efficiency after these cuts is ~ 30%. At this point, on each of the two legs there is still an 
ambiguity due to the fact that each lepton can be either the product of the first or of the second step in the 
decay chain. This gives 4 possible combinations. Furthermore, the full reconstruction results in a quartic 
equation which can have zero, two or four solutions. We show in Fig. © the distribution of the calculated 
A mass for all of the retained combinations as a full line. The dashed line shows the combinations with 
the wrong lepton assignment. A clear and narrow peak emerges over the combinatorial background. The 
width is approximately 6 GeV, determined by the resolution of the measurement of the momentum of the 
leptons. 

5. CONCLUSIONS 

In this contribution, we have described a novel technique for reconstructing the mass and momenta of 
SUSY particles. This technique does not rely on any approximate formulae nor on endpoint measure- 
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Fig. 5: Reconstructed mass of the A. 

merits. All events contribute to the sparticle mass determination and decay kinematics reconstruction, 
even if they are away from the endpoint of the distribution. The method may be particularly useful when 
the SUSY mass scale is large. In that case the statistics can be so low that the endpoint cannot be seen 
clearly while the SUSY sample itself is very clean. 

The method applies most effectively when we know some of the sparticles' masses exactly, be- 
cause the number of unknown parameters in e.g. Eq.© is reduced. In the particular case where some of 
the sparticle masses are measured at a LC the sparticle cascades may be solved completely and study of 
the decay distributions and higher mass determination becomes possible at LHC. 

When all the sparticle masses are known the neutralino momentum can be reconstructed if four 
sparticles are involved in the cascade decay. The sparticles would be pair produced, and if we can 
identify both of the cascade decay chains in the events then we only need six sparticles in the cascade 
decay to solve both of the neutralino momenta on account of the missing momentum constraint. The 
reconstruction of sparticle momenta provides us with an interesting possibility for studying the decay 
distribution at the LHC. 

On the other hand, our method is not valid when some of the particles in the cascade produce hard 
neutrinos. This is unfortunately the case when the chargino decays into (s)leptons, when a f is involved 
in the decay, or when a Wis produced and decays leptonically. If such SUSY decay processes dominate 
then this method may not be useful. 

Acknowledgments 

This work was performed in the framework of the workshop: Les Houches 2003: Physics at TeV Scale 
Colliders. We wish to thank the staff and organisers for all their hard work before, during and after the 
workshop. We thank members of the ATLAS Collaboration for helpful discussions. We have made use 
of ATLAS physics analysis and simulation tools which are the result of collaboration-wide efforts. DRT 
wishes to acknowledge PPARC and the University of Sheffield for support. 

References 

[1] ATLAS Collaboration, . CERN-LHCC-99-14/15. 

[2] I. Hinchliffe, F. E. Paige, M. D. Shapiro, J. Soderqvist, and W. Yao, Phys. Rev. D55 (1997) 
5520-5540, |hep-ph/ 9610544|. 



1 1 i 1 1 1 1 i 1 _ 

I I All combinations 




[3] I. Hinchliffe and F. E. Paige, Phys. Rev. D61 (2000) 095011, |hep-ph/ 990751 9| . 

[4] H. Bachacou, I. Hinchliffe, and F. E. Paige, Phys. Rev. D62 (2000) 015009, |hep-ph/9907518| . 

[5] B. C. Allanach, C. G. Lester, M. A. Parker, and B. R. Webber, JHEP 09 (2000) 004, 

|hep-ph/ 0070 091 . 

[6] S. Abdullin et. al, CMS Collaboration, J. Phys. G28 (2002) 469, | |hep-ph/9806366| . 
[7] M. M. Nojiri, G. Polesello, and D. R. Tovey, Contribution to these proceedings (2004). 
[8] B. C. Allanach et. al, Eur. Phys. J. C25 (2002) 113-123, |hep-ph/0202233| . 
[9] G. Corcella et. al, JHEP 01 (2001) 010, jlhep-ph/0 0113 63| . 

[10] S. Moretti, K. Odagiri, P. Richardson, M. H. Seymour, and B. R. Webber, JHEP 04 (2002) 028, 

|hep-ph/ 02 41231 . 

[1 1] E. Richter-Was, D. Froidevaux, and L. Poggioli, ATLAS NOTE ATL-PHYS-98-131. 

[12] B. Gjelsten, J. Hisano, K. Kawagoe, E. Lytken, D. Miller, M. Nojiri, P. Oslandi, and G. Polesello, 
see http://www. ippp. dur.ac. uk/^georg/lhclc. 

[13] T. A. Collaboration, CERN/LHCC/99- 14 ATLAS TDR 14 (1999). 

[14] S. Abdullin, D. Denegri, and F. Moortgat, CMS NOTE-2001/042 — PHYSICS. 



